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Scaling  Underwater  Exploding  Wires 

J.  R.  McGrath 

Energy  Conversion  Branch 
Electronics  Division 

Five-mil  Nichrome  wires  have  been  exploded  under  water  using  a  l/2-/cf  capacitor 
which  stored  energy  up  to  100  joules.  The  results  indicate  that  peak  pressure  scales  like  a 
chemical  explosive  if  losses  due  to  the  circuit  and  joule  heating  of  the  wire  itself  are 
accounted  for.  On  this  basis,  the  equivalent  weight  of  TNT  represented  by  the  electrical 
energy  stored  in  the  capacitor  is 

W  =  525(0.325CFo*  -  £,)///„ 

where  W  is  the  weight  in  micropounds,  C  is  the  capacitance  in  farads,  is  the  initial 
charging  voltage,  Et  is  the  estimated  energy  dissipated  by  the  wire  in  joules,  and  Hn  is 
the  heat  of  detonation  of  TNT  in  calories  per  gram.  This  scaling  behavior  extends  the  law 
of  similarity  six  decades  in  terms  of  weight,  from  pounds  to  micropounds.  The  peak  pres¬ 
sure  for  exploding-wire  phenomena  has  been  obtained  from  data  and  is  empirically 
expressed  as 

p„  =  26,800 


where  pm  is  peak  pressure  and  R  is  pressure-gage  distance. 

The  instrument  response  to  the  short-duration  shock  wave  gives  rise  to  a  new  defini¬ 
tion  of  the  explosion  constant  r,.  The  reduced  time-constant  parameter  shows  qualitative 
agreement  in  value  and  slope  with  chemical  data  scales,  and  is  given  empirically  as: 


ffniaya.n 


INTRODUCTION 

This  study  relates  the  research  areas  of  exploding-wire  phenomena  (EWP)  and  chemical 
underwater  explosions  (CUE).  In  this  report  the  relevant  theories  and  experiments  of  each 
held  are  outlined;  and  upon  this  relationship  the  purpose  of  this  study  is  established. 

Exploding-Wire  Phenomena 

Since  Nairn?  (1)  first  reported  his  EWP  experiments  in  1774,  considerable  research  has 
been  conducted  on  this  subject.  At  the  present  time,  most  of  the  phenomenological  aspects  have 
been  explored  in  some  depth.  In  the  last  two  decades,  considerable  information  has  been 
amassed  with  the  aid  of  improved  instrumentation,  most  notably  the  streak  and  framing  camera, 
the  Park  current  shunt,  and  the  Mach-Zehnder  interferometer.  However,  no  universally  ac¬ 
cepted  theory,  explaining  satisfactorily  the  complex  events  associated  with  exploding  wires, 
has  been  established. 

In  the  absence  of  a  theory  for  EWP,  the  discussion  of  experimental  work  must,  injustice, 
include  a  complete  list  of  investigations,  which  exceeds  the  scope  of  this  introduction.  However, 
a  comprehensive  annotated  bibliography  of  EWP  experiments,  equipment,  first-order  theories, 
and  techniques  is  given  by  Chace  and  Watson  (2).  A  historical  review  of  EWP  research  trends  has 
been  compiled  by  McGrath  (3).  A  monograph  on  exploding  wires,  drawn  from  the  papers 
presented  at  three  conferences  devoted  to  the  subject,  has  been  edited  by  Chace  and  Moore  (4). 
General  discussions  of  the  processes  taking  place  in  an  exploding  wire  by  Bennett  (5)  and 

NRI.  Problem  KOI -01 ;  ONR  Problem  RR  010-04-41-5950;  RuShips  Problem  SR  007-12-01-0800.  This  ts  an  interim  report  on  one 
phase  of  the  problem;  work  on  this  and  other  phases  is  continuing.  Manuscript  submitted  April  4,  1965. 
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Chace  (6)  give  some  idea  of  the  complexity  of  the  phenomena  and  the  increasingly  sophisticated 
techniques  used  to  observe  it. 

Review  of  Underwater  Exploaion  Theory 

In  this  ?.  <n  the  concepts,  theories,  and  experiments  of  chemical  underwater  explosions 
(CUE)  which  bear  directly  upon  this  study  are  outlined. 

For  comparison  of  one  CUE  event  with  another,  resort  is  made  in  common  practice  to 
similarity  curves.  These  curves,  in  turn,  are  based  upon  the  principle  of  similarity  which  was 
first  credited  to  Hilliar  (7).  A  more  comprehensive  discussion  of  similarity  and  scaling  is  given 
by  Snay  (8).  Essentially,  this  principle  states 

Pm{R,t)  =  Pm(nR,  nt)  (1) 

where  pm  is  the  peak  pressure,  R  is  the  pressure-gage  distance,  t  is  the  time,  and  n  is  a  constant. 
Physically,  this  means  that  if  an  explosive  charge  in  the  form  of  a  cube  having  a  side  of  length 
L  is  enlarged  to  a  side  of  length  nL,  then  the  same  pressure  will  be  measured  a  distance  nR  away. 
The  duration  of  the  explosion  wave  increases  by  the  factor  n.  The  utility  of  this  principle  is  that 
it  facilitates  the  calculation  of  shock-wave  parameters  (peak  pressure  and  time)  for  charges  of 
different  weights  (equivalendy,  different  volumes).  Theoretically,  the  principle  of  similarity 
satisfies  the  basic  equadons  of  modon  of  a  fluid  and,  for  a  shock  wave,  the  Rankine-Hugoniot 
relations.  Experimentally  the  principle  of  similarity  is  important  because  pressure  and  ume 
functions  can  be  determined  by  measurements  by  varying  either  R  or  the  charge  size  (volume). 
The  other  variable  is  determined  by  the  principle  of  similarity.  Since  the  peak  pressure  of  a 
shock  wave  decreases  faster  than  R~\  a  suitable  theory,  sadsfying  the  principle  of  similarity, 
and  accommodadng  this  experimental  observation,  is  desired. 

An  account  of  underwater  research  prior  to  World  War  II  is  given  by  Kennard  (9).  The 
advent  of  World  War  II  stimulated  American  and  British  invesdgators  (10).  Shortly  after  World 
War  II,  the  theoretical  and  experimental  work  in  this  field  was  drawn  together  by  Cole  (11). 
Several  theories  were  introduced  to  explain  shock-wave  propagation  in  water:  Penney  (1940), 
Penney-Dasgupta  (1942),  Kirkwood-Bethe  (1942),  Kirkwood-Brinkley  (1944),  and  Osborne- 
Taylor  (1944). 

In  this  study  the  asymptotic  form  of  the  Kirkwood-Bethe  theory  (12)  is  used  because  most 
work  is  compared  to  it.  This  theory  in  general  treats:  (a)  the  specification  of  the  initial  conditions 
at  the  water-explosion  gas  interface,  (b)  the  theory  of  propagation  of  a  spherical  shock  wave  in 
water,  and  (c)  the  motion  of  the  subsequent  gas  bubble.  The  asymptotic  form  of  the  Kirkwood- 
Bethe  thec.iv  leads  to  an  expression  for  peak  pressure  of  the  form: 


where  a0  is  the  explosive-charge  radius. 

In  general  this  theory  is  based  upon  the  formalism  of  the  fundamental  equations  of  hydro¬ 
dynamics.  It  does  explicitly  satisfy  the  principle  of  similarity.  Finally,  it  can  be  applied  easily 
to  a  wide  variety  of  e:  p'osives  (a  consideration  which  affects  only  the  initial  conditions). 

Purpose  of  Report 

A  review  of  the  literature  reveals  no  studies  that  have  been  made  to  determine  if  any 
similarity  exists  between  CUF.  and  EWP  events  The  purpose  of  this  report,  therefore,  is  to 
determine  if  such  similarity  exists  ind  what  the  scaling  factors  are. 
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DISCUSSION 

In  this  section  the  techniques  by  which  the  data  are  interpreted  are  discussed,  and  the  basic 
assumptions  used  throughout  the  study  are  specified.  These  specific  considerations  are:  (a)  a  cor¬ 
rection  factor  for  the  instrument  response  to  the  short-duration  shock  wave,  (b)  the  energy 
dissipated  at  the  spark-gap  switch,  (c)  energy  required  for  complete  vaporization  of  the  wire, 
and  (d)  the  electrical  energy  expressed  in  terms  of  the  equivalent  TNT  weight.  These  four  items 
are  considered  in  sequence  below. 


Duhamel’s  Theorem 


Usually,  for  large-amplitude  chemical  explosions,  the  response  of  the  detecting  instrument 
is  of  nominal  interest,  because  the  explosion-time  parameter  r,  is  large  compared  to  the  response 
time  of  the  instrument  t0  (Eqs.  4  and  5).  For  shock  waves  of  small  duration,  however,  the  varia¬ 
tion  of  pressure  with  time  can  be  appreciable  when  compared  to  the  response  time  of  the  de¬ 
tection  system.  The  result  of  this  situation  is  an  “apparent”  but  untrue  pressure-time  profile. 
To  accommodate  this  nonideal  situation,  a  mathematical  technique  is  utilized  to  correct  for  the 
physical  shortcomings  of  the  detection  system.  This  technique  is  known  as  Duhamel’s  theorem 
(13),  or  the  superposition  theorem. 

Generally  the  performance  of  a  system  may  be  detei  mined  not  only  by  its  response  to  a 
unit  step  function,  but  also  by  its  response  to  an  arbitrary  function.  These  are  related  (in  the 
notation  of  this  report)  by  the  theorem  of  Duhamel: 


A- o 


h(t)  —  p(0)A(t)  +  I  p'(k)A(t-k)d\ 


(3) 


where  h  ( t )  is  the  oscilloscope  output  response  in  volts  to  an  arbitrary  function  p  ( t ) ,  the  incident 
pressure  wave.  Here  p  (0)  is  p(t)  at  t  —  0,  k  is  a  dummy  variable,  A  (f)  is  the  response  of  the 
oscilloscope  plus  associated  components  in  volts  to  a  unit  step  function  of  pressure,  and  h(t)  is  an 
experimentally  measured  quantity  ( i.e .,  the  trace  on  the  crt  in  volts  as  a  function  of  time).  The 
character  of  p(t)  and  A  (t)  depend  upon  the  amount  and  type  of  explosive  and  on  the  details  of 
the  detection  system,  respectively. 

Equation  (3)  generally  indicates  that  if  A(t)  is  expressed  analytically  in  terms  of  parameters 
such  as  electric  current  or  mechanical  “time  constants,”  and  if  h{t)  is  also  approximated  by  some 
parameters,  then  p(t)  can  be  solved  analytically  by  using  Eq.  (3)  in  principle  by  means  of  the 
Laplace  transform.  On  the  other  hand,  if  A(t)  and  h(t)  are  given  in  tabular  form,  then  p(t)  can 
be  solved  numerically,  as  outlined  by  Osborne  (14).  The  solution  to  Eq.  (3)  by  either  method 
yields  a  factor  which  is  used  to  correct  the  observed  first  peak  voltage  ort  the  oscillogram  to  the 
peak  pressure  in  the  shock  wave.  This  peak  pressure  is  frequently  the  quantity  of  greatest  in¬ 
terest;  therefore,  let  us  assume  for  the  sake  of  simplicity  that  the  pressure  due  to  the  shock 
wave  at  a  point  is  reasonably  described  as 


p(t) 


0 

.Pmt,h' 


for  l  ^  0 
for  l  0 


(4) 


and  that  /4(t),  the  total  instrument  response,  is  approximated  by 


A(t 


-  *  •".) 


for  t  <  0 
for  t  »  0. 


(5) 


Substitution  into  and  integration  of  Eq.  (3)  using  these  definitions  yields  the  expression 


h(t) 


PmA. 


(6) 
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Of  course,  pm  is  the  peak  pressure  of  the  shock  wave  in  pounds  per  square  inch  and  Ax  is  the 
maximum  static  (t  =  *)  constant  in  volts  per  pound  per  square  inch  of  the  oscilloscope,  in¬ 
cluding  the  pressure  gage  and  amplifier  system.  The  time  at  which  h(t)  becomes  a  maximum 
value  is 


t mar 


(7) 


The  data  of  this  report  indicate  that  rf  and  ra  are  of  the  same  order  of  magnitude.  For  this  case, 
if  ye  lies  between  1/2  {r0)  and  3/2(r0) ,  Eq.  (7)  can  be  approximated  by  the  expression 


tmox  2 


(8) 


Thus,  when  the  response  time  of  the  detection  system  (oscilloscope,  amplifier,  and  pressure 
gage)  to  a  step  function  of  pressure  is  of  the  same  order  of  magnitude  as  the  decay  time  of  the 
explosion  wave,  the  observed  time  of  rise  to  peak  voltage  on  the  crt  is  approximately  equal  to 
the  decay  time  of  the  actual  (or  hypothesized)  explosion  wave.  The  value  of  hit)  at  its  maximum 
point  is  given  by  substituting  Eq.  (7)  into  Eq.  (6),  which  yields 


h(t)max  —  PmAx  {(t  j  [exp“(V<r«'TaW|n<r^Ta>— exp-(V(Te'To»ln(Te/,a)]  j.  (9) 

As  a  result,  a  measured  peak  deflection,  h{t)max  in  volts,  on  an  oscillogram  under  these  condi¬ 
tions  may  be  corrected  to  yield  the  true  peak  shock-wave  pressure  in  pounds  per  square  inch  by 


Pm  = 


hit). 


4X  ^ ■-  j  [exp-,Ta'Tf-Ta>ln<Vf«>  -  exp-(VTe-fa>ln<VTa>]  j 


(10) 


If  pm  is  defined  as  A(t)Mar/4x,  then  Eq.  (10)  may  be  written 


Pm  ~ 


Pm 


[exp-(Ta,V-ra>ln<Vra>  -  exp“(rdV-Ta»n<Tdra>] 


(ID 


I  o  the  extent  to  which  Eqs.  (4)  and  (5)  are  reasonable  approximations,  so  far  as  peak-pressure 
measurements  are  concerned,  the  comparatively  slow  response  of  the  detection  system  to  a 
rapidly  varying,  short-duration  pressure  wave  is  overcome.  Equations  (8)  and  ('0)  give  the 
principal  explosion-wave  parameters  of  interest,  p„  and  r,.  directly  from  the  observed  oscillo¬ 
gram  trace,  hit).  The  quantity  pZ  may  be  thought  of  as  a  peak  pressure  taken  from  the  maxi¬ 
mum  oscillogram  deflection,  uncorrected  for  the  finite  response  time  of  the  apparatus. 


Circuit'  Switch  Loss 

According  to  Moses  and  Kornefl  in  a  recent  paper  (4a),  an  air-gap  switch  produces  an 
initial  current  slope  of  zero  and  inhibits  the  current  rise  in  the  wire  to  a  value  less  than  that 
predicted  by  circuit  theory  during  the  first  quarter  cycle.  Their  study  shows  that  approximately 
35  percent  of  the  initial  stored  energy  of  the  capacitor  is  dissipated  by  the  switch  during  the 
first  ha'f  he  disc  harge.  Electrical  energies,  assuming  no  circuit  loss,  are  designated  E, 

w  here  i  jl  energy  of  the  capacitor;  electrical  energies  accounting  for  circuit  losses  and 

subsequent  assumptions  are  designated 
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Energy  Required  to  Heat  Wire 

Energy  is  required  (o  vaporize  completely  (an  assumption)  the  wire.  This  energy  require¬ 
ment  is  taken  into  account  by  Vaughan  (15)  using  the  handbook  values  of  the  thermodynamic 
properties  of  the  element  by  the  general  expression 


Ec  —  kLD*  (12) 

where  Ev  is  the  vaporization  energy  in  joules,  k  is  a  material-dependent  constant,  L  is  the  length 
of  the  wire  in  inches,  and  D  is  the  wire  diameter  in  mils.  For  Nichrome  wire  (used  throughout 
this  study),  k  =  0.572  in.'1  mils'1  joules1.  The  energy  assumed  required  for  complete  vaporiza¬ 
tion  of  half-inch-long,  five-mil-diameter  Nichrome  wire  is  estimated  from  Eq.  (12)  to  be  a 
constant:  ten  joules. 

Weight  Equivalence  of  Electrical  Energy 

The  electrical  energy  initially  stored  in  the  capacitor  may  be  expressed  in  terms  of  an 
equivalent  weight  of  TNT  by  means  of  the  heat  of  detonation.  For  TNT,  Price  (16)  indicates  the 
calorimetric  heat  of  detonation  to  be  1013  calories  per  gram.  Consequently,  the  equivalent 
weight  W  in  micropounds  is 


rrtr  / 0.325 CVl-EA 

r_525  l — m, — )  ('« 

where  C  is  the  capacitance  in  farads,  V0  is  the  initial  voltage  in  volts  to  which  the  capacitor  is 
charged,  £‘r  is  estimated  energy  in  joules  required  to  vaporize  the  wire  completely,  and  Hn  is 
the  heat  of  detonation  of  TNT  in  calories  per  gram.  The  equivalent  weight  W  also  includes  the 
circuit  loss  of  35  percent  at  the  spark-gap  switch. 

EXPERIMENTAL  ARRANGEMENT 
Equipment 

The  equipment  employed  in  this  study  is  divided,  according  to  its  use,  into  three  general 
systems:  (a)  a  charging  circuit,  (b)  a  discharge  circuit,  and  (c)  a  detection  system.  A  discussion 
of  the  operation  of  this  equipment  follows  a  brief  description  of  each  system,  and  its  arrange 
tnent  is  shown  in  Fig.  1. 


Lo 

-dWSRT — ^ 


POWER  *- 
S  JPPLY 


GAGE- 

GABlr 

UNIT 


AMPLIEIER 


lO 


OSCILLOSCOPE 


TO  VACUUM  PUMP 
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*■  «X-  s  -  '  *>*•  |»‘Hrt  supple  (Ml)  jik!  voltmctri  (right) 


I  he  series  charging  . . •  <  t< it  re. twists  of  the  following  c  nni{>onents:  (a)  a  (lower  sit j>ji!\,  {!>)  a 
charging  resistor,  (c)  a  low -indue  tatue  capacitor,  and  (d)  a  soheneter  plat  ed  in  parallel  with  the 
capac  itor.  The  energy  stored  m  the  capacitor  is  determined  hv  the  voltage  to  which  it  is  charged 
bv  the  power  supply  (Figs.  2.  3). 

1  he  discharge  circuit  is  a  series  1,R(.  circuit  consisting  oi  the  following  toniponcnts  (a)  the 
low -inductance  capacitor,  (b)  a  sficcial  air-gap  swttth  (17).  (t)  elettn«les,  (d)  vacuum  pump, 
and  (e)  four  low-inductance  k,.cls  ft om  the  spark-gap  switch  to  the  electrodes  (Fig.  4).  The 
dielectric  strength  of  die  electrical  cable  is  sufficient  to  withstand  fifVHMI  soils. 
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Fig.  4  —  The  discharge  circuit  (vacuum  pump  below  capacitor,  not  shown) 


The  detection  system  consists  of  the  following  components:  (a)  an  oscilloscope,  (b)  an 
oscilloscope  camera,  (c)  an  oscilloscope  plug-in  amplifier,  and  (d)  a  tourmaline  pressure  gage. 

Methods 

Operation  of  the  discharge  circuit  is  restricted  to  20,000  volts,  since  breakdown  at  the  air- 
gap  switch  occurs  at  23,000  volts  at  atmospheric  pressure.  When  the  capacitor  is  charged  to 
some  voltage  (To),  the  electrical  energy  stored  in  the  capacitor  may  be  switched  through  the 
LRC  circuit  by  evacuating  the  spark-gap  switch  with  the  vacuum  pump.  The  current  passes 
through  the  low-inductance  leads  to  the  electrodes  in  the  water  tank,  and  back  to  the  capacitor. 

The  capacitance  of  the  low -inductance  capacitor  was  measured  by  an  impedance  bridge 
to  be  0.507  ±  0.001  p.f.  The  internal  inductance  was  stated  by  the  manufacturer  to  be  60  milli- 
microhenries.  The  ringing  frequency  of  the  circuit  was  observed  to  be  about  400  kc,  and  the 
circuit  inductance  was  0.31  microhenry .  The  circuit  resistance  was  0.05  ohm.  Evaluation  of 
circuit  parameters  required  the  use  of  oscillograms,  and  the  analytic  technique  was  adopted 
from  Page  (18).  The  time  required  for  the  current  oscillations  to  damp  out  was  about  25  /isec 
when  the  test  gap  was  shorted  with  a  copper  bar. 

The  sweep  speed  of  the  oscilloscope  was  kept  constant  at  five  microseconds  per  centimeter 
throughout  the  study.  The  time  count  was  initiated  at  the  oscilloscope  by  placing  a  pickup  loop 
near  the  spark-gap  switch.  When  the  discharge  began,  a  voltage  was  induced  in  the  loop, 
and  scope  time  was  counted  from  that  instant.  Because  the  transit  *ime  of  the  shock  wave  from 
the  exploding  wire  to  the  pressure  gage  varied  tine  to  the  variation  of  the  gage  distance,  a 
time  delav  was  employed  on  the  oscilloscope  in  order  to  present  the  interval  in  which  the  shock 
wave  passed  b\  the  gage.  The  value  of  the  time  delav ,  of  course,  depended  upon  the  distance  be¬ 
tween  the  gage  and  the  wire.  The  sweep  time  was  c  alibrated  In  means  of  a  time-mark  generator. 

The  presentation  of  the  pressure-time  history  of  the  shock  wave  on  the  ert  of  the  oscillo- 
scojte  was  recorded  by  the  Polaroid  oscilloscope  camera.  Each  oscillogram  was  then  enlarged, 
and  a  tracing  of  each  enlargement  was  made.  Emm  these-  enlargements,  measurements  of 
shoc  k-wave  parameters  (amplitude  and  time)  were  made. 

I  he  differential  .ouplifier  was  used  to  exclude  some  extraneous  electromagnetic  signals. 
The  rise  time  of  the  oscilloscope  and  the  amplihoi  as  i  unit  was  about  0.2  (manufac¬ 
turer’s  data).  The  amplihet  was  calibrated  lor  each  se'ting  used  (vertical  vale  0.5  volt  cm. 
I  volt  cm,  and  2  volts  cm  hv  using  foui  or  live  known  signals).  I  bc  amplihet  response  was 
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in  error  by  a  linear  factor.  The  measured  pressures  in  Eq.  (15)  were  corrected  for  the  ampli¬ 
fier’s  error. 

The  pressure  gage  used  in  this  study  was  a  commercially  available,*  l/4-in.-diameter, 
4-ply,  underwater  type,  tourmaline  pressure  gage.  The  design,  use,  and  performance  of  such 
gages  have  been  studied  by  Arons  and  Cole  (19),  and  also  by  Fox,  et  al.  (20).  The  unit  was 
mounted  on  16  ft  of  shielded  cable  which  ter  minated  directly  at  the  amplifier  inout.  The  man¬ 
ufacturer’;.  crystal  constant  (hydrostatically  determined)  was  1.99X10-12  c  .omb/psi.  The 
gage  was  employed  in  an  edge-on  position  to  the  exploding  wire.  The  capacitance  of  the  gage- 
cable  unit  was  423.7  ±  0.1  fifif.  The  impedance  for  the  feed-cable  length  was  69  ohms.  Based 
upon  these  figures,  the  estimated  delay  time  for  the  feed  cable  was  about  29  nanoseconds,  not 
a  source  of  concern  in  this  experiment. 

To  insure  that  the  gage  was  watertight,  a  thin  coat  of  C-276f  was  applied.  X-ray  photo¬ 
graphs  of  the  gage  in  the  modified  configuration  indicate  that  the  coating  was  0.5  mm  thick 
in  front  of  the  edge-on  position.  Moreover,  the  diameter  of  the  gage  was  measured  to  be  7.65 
mm. 

To  insure  that  the  gage-cable  unit  was  functioning  properly  after  the  exploding-wire  ex¬ 
periment,  calibration  and  field  tests  were  performed  at  the  Naval  Propellant  Plant  by  Slifko.t 
The  gage-cable  unit  was  calibrated  by  the  bursting-diaphragm  method  which,  in  the  modified 
form,  yielded  a  crystal  constant  of  2.13  X  10~12  coulombs/psi,  checking  to  within  7  percent  of 
the  manufacturer’s  value.  In  additi:  n,  tests  of  the  gage  in  the  edge-on  position  were  conducted 
using  standard  Army  Engineer  detonator  caps  (about  a  millipound  of  TNT).  These  tests  indi¬ 
cated  that  the  gage  performed  in  a  satisfactory  manner,  one  test  is  compared  to  an  EWP  event 
in  Fig.  5.  Actually,  the  rise  time  for  the  gage  was  a  microsecond  faster  than  other  gages  against 
which  it  was  compared  ( i.e about  four  as  compared  to  five  microseconds).  The  pressure  mea¬ 
surements  were  corrected  to  reflect  the  gage  constant. 


Fig.  5  —  Comparison  of  a  CUF.  and  underwater  FiWP  event  using 
normalized  pressure  irrmi  lime.  CUF,  data  courtesy  of  John  P.  Slif  *  i. 
NOL  Unit,  Indian  Head,  Md. 


•Crystal  Research,  Inc..  Cambridge,  Mass. 
tZophar  Mills,  Brooklyn,  N.t . 

tjohn  P.  Slilko,  NOl.  Unit.  Naval  Propellant  Plant,  Indian  Head.  Mil 
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Based  upon  the  previous  discussion,  the  voltage  (vertical  deflection)  measured  from  the 
oscillogram  is  given  by  Cole  (1  la)  as 


KA 

Co  +  Ci 


) 


(14) 


where  V  is  the  yoltage  in  volts,  KA  is  the  crystal  constant  in  micromicrocoulombs/psi,  Co  is  the 
lumped  capacitance  of  the  cable-gage  unit  in  micromicrofarads,  C 1  is  the  input  capacitance  of 
the  amplifier  in  micromicrofarads,  and  p*  is  the  peak  pressure  in  psi  uncorrected  by  Duhamel’s 
theorem.  Since  the  only  quantity  unc.  termined  is  p*,  the  value  of  p*  is,  for  Co  =  424  ppf, 
Ci  =  45  ppf,  and  KA  —  2.13  /A/xc/psi: 


p*  (psi)  =  220  v.  (15) 

Equation  (15)  does  not  reflect  the  correction  required  for  the  linear  amplifier  error  nor  due  to 
the  system  response  (Eq.  10). 

The  distance  between  the  pressure  gage  and  the  exploding  wire  was  varied  from  1.5  to 
12.0  in.;  measurements  were  taken  every  1.5  in.  Because  the  distances  involved  were  small, 
and  because  accurately  positioning  the  gage  was  not  possible,  accurate  distance  measurements 
were  difficult  to  make.  The  actual  distance  was  always  less  than  the  measured  distance  by  as 
much  as  0.200  in.  due  to  the  technique  of  mounting  the  wire  and  measuring  the  distance. 

DATA  AND  RESULTS 

The  data  used  in  this  study  are  obtained  f  rom  oscillograms  which  give,  directly,  the  voltage 
deflection  from  a  reference  position  on  a  crt  as  a  function  of  time.  The  voltage  deflection  results 
from  the  charge  applied  to  the  plates  of  the  crt.  This  charge  is  generated  by  the  piezoelectric 
effect  plus  the  amplifier;  in  this  case,  it  is  due  to  the  compression  of  a  tourmaline  gage  by  a 
shock  wave.  The  data  in  this  study  are  obtained  from  the  use  of  Nichrome  wire  exclusively. 

A  comparison  of  EWP  pressure-time  profiles  can  be  made  to  CUE.  In  Fig.  5,  for  example, 
a  profile  representing  a  CUE  and  one  an  underwater  EWP  event  a*  e  presented  to  show,  qualita¬ 
tively,  the  contrast  between  the  two  types  of  explosions.  The  pressure  of  each  is  normalized  for 
this  comparison. 

As  mentioned  before,  and  throughout  this  article  for  EWP  data,  a  distinction  is  maintained 
between  the  peak  pressure  pni,  which  is  measured  directly  from  oscillographic  data,  and  p,„, 
which  is  calculated  using  p*  on  the  basis  of  certain  time  constants.  The  quantity  p *  includes  all 
corrections,  calibrations,  and  assumptions  except  those  given  by  the  Duhamel  theorem.  The 
quantity  p,„  includes  all  corrections,  including  those  given  bv  DuhaimTs  theorem,  i.e.,  the 
denominator  in  Eq.  (10)  or(  I !). 

Pressure  Measurements 

Energy  Adjustments  Related  to  Scaling  Peak  Pressure  p*m—  In  CUE  experiments,  the  peak 
pressure  shows  its  distance  dependence  in  the  asymptotic  form  as 

Pm  «/?-*«  (16) 

where  R  is  the  gage  distance  and  /3o  is  a  constant  slightly  greater  than  unity  for  CUE  data. 
The  peak  pressure  also  depends  upon  the  charge  weight  W.  At  any  value  of  R  the  peak  pres¬ 
sure  is  related  to  the  weight  of  the  charge  anil  the  gage  distance  R  by  the  law  of  similarity: 
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P.  -  r )  ■  <17) 

Equivalently  the  electrical  energy  of  the  capacitor  could  be  related  by  Eq.  (17)  through  the  heat 
of  detonation  Ho  to  an  expression  involving  chemical  energy.  A  similarity  curve  for  pm  versus 
either  k\  (Wl,3IR)*>  or  Jt2  {Eil3IR)%>  is  possible.  The  literature,  however,  historically  uses  the 
relation 


/jpi/3\0o 

=  (18) 

where  k i  depends  upon  the  type  of  explosive. 

In  Fig.  6  the  average  peak  pressure  p,*  in  arbitrary  units  is  plotted  against  the  cube  root  of 
the  corrected  ( ECOrr )  and  uncorrected  ( E )  initial  capacitor  energy  at  a  fixed  distance.  Clearly 
the  corrected  electrical-energy  values  form  a  straight  line  whose  slope  is  greater  than  that 
for  the  acoustic-law  curve.  The  uncorrected  energy  values  do  not  form  a  straight  line.  All  wires 
used  were  5-mil-diameter,  one-half-inch-long  Nichroine. 

Oscillographic  data  were  taken  for  capacitor  energies  varying  from  25.4  joules  (10  kvj  to 
101.5  joules  (20  kv)  and  for  gage  distances  varying  from  1.5  in.  to  12.0  in.  in  increments  of  1.5 
in.  The  results  of  this  experiment  are  shown  for  peak  pressure  in  Fig.  7.  Here  the  peak  pressure, 
p*,  is  plotted  against  the  cube  root  of  the  equivalent  weight  of  TNT  divided  by  the  gage  distance. 


ELECTRICAL  ENERGY  E^5  (  joules'7-4  ) 


Fig.  6  —  Peak  pressure  (uncorreclcd  for  instrument  res|>onse)  versus 
the  eube  root  of  the  electrical  energy  (corrected  and  uncorrccted  for 
circuit  losses)  for  Nichroine  wire,  1/2  in.  long,  5  mils  in  diameter,  at 
a  gage  distance  of  9  in. 
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The  equivalent  weight  of  TNT  is  estimated  for  100  percent  of  1/2  CVl  and  for  100  percent  of 
1/2  CV$  less  the  values  of  the  circuit  loss  (35  percent  of  1/2  CVh  and  wire  vaporization  (lOjoules). 
These  considerations  are  termed  uncorrected  ( E )  and  corrected  electrical  energies  (£Wr), 
respectively.  The  corrected  data  in  Fig.  7,  therefore,  represents  an  expression  of  the  form 


=k<hr ) 


(19) 


where  p*  here  is  the  peak  pressure  uncorrected  by  Duhamel’s  theorem,  k3  is  a  constant  of  the 
same  order  of  magnitude  as  k->  for  TNT,  k*  is  a  constant,  ECorr  is  defined  above,  and  /3 1  is  a 
constant  slightly  greater  than  unity  for  EWP  data.  Figure  7  contains  data  for  the  two  extreme 
values  of  electrical  energy  only,  25.4  and  101.5  jouies.  Note  that  the  data  using  only  corrected 
energies  form  a  common  curve,  or  scale,  whereas  the  uncorrected  energies  form  two  distinct 
lines  having  the  same  slope.  The  data  for  Ecnrr  is  chosen  to  represent  EWP  because  it  forms  a 
common  curve.  On  the  basis  of  the  information  presented  in  Fig.  7,  further  comparisons  are 
made  using  the  corrected  electrical  energy  values  in  Fig.  8. 

Comparison  of  EWP  p, to  CUE  pm  Data  —  I n  Fig.  8  CUE  and  EWP  data  are  compared.  In¬ 
dividual  CUE  data  points  arc  included  (see  legend)  from  the  work  of  several  investigators 
Two  curves  ate  presented  relating  peak  pi  ssure,  pm  for  CUE,  and  p$  for  EWP,  to  the  scaled 
reciprocal  distance.  Curve  1  is  the  asympt  aic  approximation  of  the  Kirkwood-Bet  he  theory 
using  TNT  (density  1.59  g/cc)  and  is  given  by  the  empirical  expression  reported  by  Arons  (24) 
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Fig.  8  —  Peak  pressure  vrrsus  scaled  reciprocal  distance  for  CUE  data  (curve  1) 
and  EWP  data  (curse  2).  EVVP  data  arc  uncorrected  for  instrument  response 
(Duhamel's  theorem) 


(1^1/3X1.13 

r)  (20) 

where  W  is  the  charge  weight  in  pounds,  R  is  the  gage  distance  in  feet,  and  pm  is  the  peak  pres¬ 
sure  in  pounds  per  square  inch  for  CUE  events.  This  curve  is  drawn  from  data  using  55,  25,  and 
0.5  lb  TNT  charges.  Curve  2  is  the  result  of  plotting  all  exploding-wire  data  (i.e.,  20,  18,  16, 
14,  12,  and  10  kv  shots  at  distances  of  1.5,  3.0,  4.5,  6.0,  7.5,  9.0,  10.5,  and  12.0  in.)  and  adjusting 
each  energy  for  the  assumption  of  complete  vaporization  and  for  the  energy  dissipated  at  the 
spark-gap  switch.  Clearly,  several  statements  seem  to  be  in  order  upon  examining  Fig.  8.  First, 
curve  2  (EWP  data)  is  about  50  percent  lower  than  curve  1  (TNT  data)  for  a  given  value  of 
W'I*IR.  Secondly,  on  the  basis  of  Fig.  8  the  corrected  EWP  data  seems  to  form  a  similarity  cut  ve 
not  unlike  TNT  data.  Thirdly,  the  slope  of  curve  2  seems  to  be  less  than  1.13,  as  in  TNT  data, 
but  also  is  no  less  than  unity.  Curve  2  of  Fig.  8  supports  the  conclusion  that  the  peak  pressure 
at  least  of  EWP  scales  like  a  chemical  explosion  in  general,  when  the  electrical  energy  is  cor¬ 
rected  for  circuit  and  evaporation  requirements. 

The  weight  of  an  equivalent  TNT  charge  for  EWP  is  calculated  to  be 


W  (micropound)  =525 


325  CVl-E „ 


H„ 


)■ 


(21) 


If  this  equivalent  weight  of  TNI  is  cast  into  spherical  form,  its  radius  a«  in  inches  is  given 


*■**&&?<?*  - 
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a.  (in.)  =  0.0394  [(£„„)] 1/3 

where  ECorr  is  the  corrected  electrical  energy  in  joules. 

For  voltages  less  than  10  kv  (25.4  joules),  the  basic  assumption  of  complete  vaporization 
tends  to  be  less  believable  because  there  would  be  very  little  energy  left  for  the  shock  wave  if 
a  constant  10  joules  is  used  to  vaporize  the  wire. 

The  essential  pressure  and  time  data  for  all  shots  of  this  study  are  presented  in  Table  1. 
Column  1  of  Table  1  identifies  the  particular  shot,  column  2  lists  the  voltage  in  kilovolts,  and 
column  3  lists  the  gage  distance  in  inches.  Column  4  lists  the  peak  pressure  p*  in  volts.  Column  5 
presents  the  peak  pressure  p*  in  psi,  calculated  from  Eq.  ( 1 5).  The  crystal  constant  as  determined 
by  Slifko,  KA  =  2.3  pp  coul/psi,  is  used  here.  Columns  6  and  7  are  explained  later  in  this  report. 

Comparison  of  EWP  (pm)  to  CUE  (pm)  Data— Based  upon  the  values  of  tm  and  the  value  of 
Ta,  the  maximum  value  of  pm  may  be  calculated  from  Eq.  (10).  Equation  (10)  varies  little  if  indi¬ 
vidual  values  of  tm  are  not  greatly  different.  The  result  is  nearly  a  constant,  i.e.,  pm  =  p* '0.43. 
Therefore,  the  oscillographic  pressure  values  may  be  corrected  by  multiplying  each  by  the  factor 
2.32.  These  values  are  presented  in  column  6  of  Table  1  and  are  also  compaied  to  CUE  theory 
and  data  (24)  and  to  a  one-megaton  nuclear  w  eapon  (25)  in  Fig.  9.  Column  7  of  Table  1  presents 
the  values  of  scaled  reciprocal  distance  based  upon  the  equivalent  weights  W  of  TNT  cor¬ 
responding  to  the  corrected  initial  capacitor  energy  (ECOrr)-  Figure  9  indicates  that  EWP  pressure 
data  are  approximately  40  percent  larger  than  CUE  data  for  the  same  value  of  Wlli/R  and  can 
be  described  approximately  by  the  expression: 
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pm  =  26,800 


/Jp/3\  1.08 

nr) 


(23) 


Conversion  From  Electrical  to  Shock-Wave  Energy 

The  energy  of  a  shock  wave  propagating  under  water  is  given  approximately  by  Cole  (1  lb): 


47t  R1  fT 

W(R)=^-jr\  [p(t)?dt. 
Po  Cn  JnH) 


(24) 


The  calculations  are  presented  in  Fig.  10  for  energies  ranging  from  25.4  joules  (10  kv)  to  101.5 
joules  (20  kv)  at  a  constant  distance,  R  -  1.5  in.  In  this  figure  the  ratio  of  shock-wave  energy  to 
initial  stored  electrical  energy  in  percent  is  plotted  against  the  initial  stored  capacitor  voltage. 
The  percentage  ranges  from  5 1  to  62  percent  and  has  an  average  value  of  58.0  percent.  Bennett 
(26)  using  4  and  5  mil  copper  wires,  at  energies  ranging  from  20  to  37  joules,  indicates  that 
approximately  40  to  60  percent  of  the  initially  stored  electrical  energy  is  deposited  in  the  wire 
and  is  available  for  subsequent  fluid  motion.  Buntzen  (4b),  using  5-mil,  one-inch-long  copper 
wires,  exploded  underwater,  finds  the  shock-wave  energy  about  63  percent  of  that  initially 
stored  in  the  capacitor.  The  energy -conversion  factor  determined  in  tiiis  study  shows  reasonable 
agreement  with  the  independent  results  of  these  investigators,  i.e.,  an  average  value  of  58.2 
percent  approximately. 

Time  Measurements 

Time  Constant  of  Underwater  EWP  Events—  In  general,  a  comparison  of  time  measure¬ 
ments  for  CUE  and  EWP  events  show  that  the  explosion  time  constant  is  smaller  by  a  fac  tor  ol 
103  in  general  for  EWP  In  most  CUE  experiments,  for  example,  the  decay  time  constant  is 
expressed  in  milliseconds.  In  this  study,  however,  a  measurement  of  only  a  lew  microseconds  is 
observed. 
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Table  1 

Pressure  and  Distance  Data  for 
Underwater  Exploding  Nichrome  Wires 


Voltage 

(kv) 


:  Pm  W»IR 

si)  (psi)  (lb'^/ft) 


20 

1.5 

8.70 

2297 

5329 

18 

8.20 

2165 

5023 

16 

7.36 

1943 

4508 

14 

6.42 

1694 

3930 

12 

5.32 

1404 

3257 

10 

4.10 

1100 

2552 

20 

3.0 

4.5 

1198 

2779 

18 

4.16 

1107 

2568 

16 

3.76 

1009 

2341 

14 

3.26 

875 

2030 

12 

2.76 

747 

1733 

10 

2.10 

578 

1341 

20 

4.5 

2.98 

794 

1842 

18 

2.75 

726 

1684 

16 

2.50 

660 

1531 

14 

2.20 

581 

1348 

12 

1.88 

501 

1161 

10 

1.40 

376 

872 

20 

6.0 

1.94 

517 

1199 

18 

1.97 

485 

1125 

16 

1.87 

464 

10';  7 

14 

1.68 

418 

970 

12 

1.43 

356 

826 

10 

1.12 

280 

650 

20 

7.5 

1.92 

481 

1116 

18 

1.79 

449 

1041 

16 

1.62 

406 

942 

14 

1.41 

353 

819 

*  O 

*  mm 

1.20 

304 

705 

10 

0.88 

223 

517 

20 

9.0 

1.62 

406 

942 

18 

1.48 

372 

863 

16 

1.30 

326 

756 

14 

1.16 

291 

675 

12 

1.00 

253 

587 

10 

0.75 

191 

443 

20 

10.5 

1.39 

342 

793 

18 

1.25 

314 

728 

16 

1.13 

286 

664 

14 

0.97 

245 

568 

12 

0.77 

1% 

455 

10 

0  58 

150 

348 

0.1816 

0.1544 

0.1208 


0.1028 


0.0685 

0.0005 

0.0515 

0.0403 

0.0618 

0.0568 

0.0514 

0.0454 

0.0386 

0.0302 

0.0494 

0.0454 

0.0411 

0.0363 

0.0309 

0.0242 

0.0412 

0.0379 

0.0343 

0.0303 

0.0257 

0.0201 

0.0353 

0.0325 

0.0294 

0.0259 

0.0221 

0.0173 


(Table  (  ontiinies) 
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Table  1  (Continued) 

Pressure  arid  Distance  Data  for 
Underwater  Ej;ploding  Nichrome  Wires 


Shot 

No. 

Voltage 

(kv) 

R 

(in.) 

n* 

"  m 

(volts) 

p* 

(psi) 

Pm 

(psi) 

(lb'«/ft) 

43 

20 

12. C 

1.25 

314 

729 

0.0309 

44 

18 

1.14 

289 

671 

0.0284 

45 

16 

1.04 

263 

610 

0.0257 

46 

14 

0.89 

225 

522 

0.0227 

47 

12 

0.75 

191 

443 

0.0193 

48 

10 

0.59 

152 

353 

0.0151 

Fig  '.I  -  t'auup.mstin  ol  ('l  l  (turvc  1)  and  KWP  (turve  2)  data  using  peak 
pressure  and  stak'd  irtiprotal  distance  F.WP  data  arc  turret  led  Itu  inslruineni 
i  espouse. 
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Fig.  10  —  Ratio  of  shock-wave  energy  to  initial  capacitor 
energy  versus  initial  capacitor  voltage 


Definition  of  t,  for  EWP  Events— For  large-amplitude  CUE  events,  the  pressure  behaves 
approximately  as  p(t )  =  p0  e~llre,  and  the  explosion  constant  r,  is  evaluated  over  the  range 
where  the  disturbance  has  fallen  to  (1  —  e"*)  of  its  maximum  value.  On  a  semilog  plot  of  peak 
pressure  versus  time,  re  is  the  slope  (in  milliseconds)  which  is  a  constant  over  the  range  of  pm 
to  approximately  0.67pm. 

The  small-amplitude  explosions  of  this  study  (EWP),  on  the  other  hand,  do  not  lend  them¬ 
selves  to  such  an  analysis.  The  rise  of  the  pressure  wave  to  a  maximum  value  does  not  behave 
like  a  step  function  of  pressure  which  then  decays  exponentially.  Moreover,  the  pressure  decays 
almost  as  rapidly  as  it  rises  to  a  maximum.  If  the  definition  of  the  explosion  constant  for  large- 
duration  explosions  is  used  for  the  small-duration  waves  in  this  study,  then  points  on  the  oscil¬ 
lographic  record  are  chosen  in  every  instance  in  the  situation  where  the  explosion  wave  has  not 
entirely  crossed  the  crystal  area.  The  explosion  time  constant  defined  for  EWP  and  for  this 
study  is  approximated  by  Eq.  (8)  and  is  defined  as  the  average  value  of  time  over  which  this 
approximate  expression  is  valid  (i.e.,  tmax  —  1/2(t„  4-  rr). 

Definition  and  Calculation  cf  t„  for  EWP  Events— The  quantity  r„,  as  expressed  in  Eq.  (5), 
for  the  amplifier  alone,  excluding  the  gage,  is  about  0.2  psec.  A  further  consideration  is  that 
the  crt  cannot  display  a  voltage  deflection  from  its  reference  position  until  the  tourmaline  gage 
has  been  sufficiently  compressed.  Consequently,  the  transit  time  required  for  the  explosion 
wave  to  compress  the  gage  sufficiently  must  be  taken  into  account.  The  transit  time  is  defined 
here  to  be  the  time  required  for  the  explosion  wave  to  traverse  (1  —  e~')  of  the  diameter  normal 
to  the  explosion  wave.  From  x-ray  photographs,  d  is  measured  to  lie  7.65  mm.  The  time  interval 
for  the  explosion  wave  to  traverse  the  distance  (1  —  e"1)  d  is  estimated  to  be:  (0.765)  (0.63)/ 
0.15  X  10®  =  3.21  psec.  The  combined  characteristic  time  constant  for  the  gage-amplifier- 
oscilloscope  system  is  calculated  to  be,  therefore,  3.21  -I-  0.20  =  3.41  /lcscc  =  t„. 

Time  Measurements  of  EWP  Shock  Wave  Profiles—  Measu cements  of  the  partial  rise  lime 
0i,  as  shown  in  Fig.  1 1,  and  of  subsequent  time  intervals  corresponding  to  p  =  0.5  pmnx  (/>., 
0-t,  and  6s  and  to  p  —  0.25  pmax  (i-e.,  6A  and  05)  have  been  taken.  The  average  value  for  all  <?i 
measurements  is  2.07  //.sec,  for  tmnx  is  3.42  /xwc,  for  d<  is  3.12  ju.sec,  for  03  is  5.08  //sec,  for  0j  is 
4.15  /xsec,  and  for  05  is  6.92  /zsec.  These  values  are  obtained  from  Table  2,  which  gives  all  time 
measurements. 

Comparison  of  (T,)raic  With  (rr)ex^—  As  stated  before,  T,.  =  2tmax  —  T„.  The  calculated  value  of  rn 
is  3.41  psec.  The  measured  average  value  of  tmnx  is  3.42  psec.  This  measure  of  agreement  gives 
some  confidence  to  the  assumption  that  Tr=2rmnx—  r„  and  to  the  consequences  of  Eq.  (10)  or  (1 1). 

Increasing  Value  of  r.  With  Propagation—  Arons  (24)  has  shown  evidence  that  r,.  for  large- 
amplitude  shock  waves  increases  in  value  with  propagation  distance  for  charge  weights  varying 
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Fig.  1 1  —  The  form  of  an  ideal  oscillographic 
record,  indicating  the  various  time  measurements 


Table  2 


Time  Measurements  from  Oscillographic  Records 


Shot 

No. 

Voltage 

<kv) 

R 

(in.) 

0. 

(Msec) 

0, 

(Msec) 

03 

(Msec) 

0, 

(Msec) 

05 

(Msec) 

tmax 

(Msec) 

20 

12.0 

2A\ 

2.72 

5.17 

3.79 

7.33 

4.57 

10.5 

2.50 

2.76 

5.52 

3.88 

7.41 

4.22 

9.0 

2.33 

2.85 

5.43 

3.88 

7.24 

-- 

25 

7.5 

2.41 

2.76 

5.43 

3.79 

7.24 

4.22 

19 

6.0 

2.38 

2.93 

5.52 

7.33 

4.05 

13 

4.5 

2.07 

3.02 

5.17 

4.14 

6.98 

345 

H 

3.0 

2.16 

3.10 

5.00 

4.22 

6.72 

3.19 

H 

1.5 

1.72 

3.45 

5.26 

4.40 

6.72 

2.93 

B 

18 

12.0 

2.30 

2.78 

5.38 

m 

.3.97 

38 

10.5 

2.21 

2.98 

5.47 

4.03 

■El 

3.79 

32 

9.0 

1.92 

2.78 

5.57 

3.84 

7.30 

3.79 

20 

6.0 

2.02 

3.46 

5.18 

4.22 

6.72 

2.93 

14 

4.5 

2.0- 

3.36 

4.99 

4.42 

6.91 

2.93 

8 

3.0 

2.11 

3.07 

4.99 

4.13 

6.91 

3.19 

2 

1.5 

1.82 

3.26 

4.80 

4.42 

6.82 

2.59 

45 

!6 

12.0 

1.99 

274 

5.58 

3.69 

7.28 

4.33 

39 

10.5 

2.17 

2.55 

5.10 

3.59 

6.99 

4.23 

33 

9.0 

2.08 

2.65 

5.29 

3.69 

7.18 

4  13 

27 

7.5 

1.99 

3.59 

5.20 

4.63 

7.18 

3.46 

21 

6.0 

1.99 

3.40 

5.10 

4.44 

6.90 

3.17 

15 

4.5 

2.08 

3.40 

5.10 

4  35 

6.99 

3.46 

9 

3.0 

2.08 

2.27 

4.82 

4.16 

6.62 

3.26 

3 

1.5 

1.89 

3.12 

4.73 

4.25 

6.62 

3  08 

40 

14 

12.0 

2.2! 

3  46 

5.38 

4.42 

7.20 

40 

10.5 

2.21 

3.36 

5.18 

4.32 

6.91 

— 

34 

9.0 

2.11 

3.36 

5.09 

4.22 

6  91 

3.46 

28 

7.5 

2.02 

3  46 

5  38 

4.42 

7.20 

3.46 

22 

6.0 

1  82 

3.36 

5.09 

4.42 

6.91 

3.08 

16 

4.5 

2.02 

3.17 

4.80 

4.32 

6  82 

10 

3.0 

1.92 

3.36 

5.09 

4.32 

6.82 

3  17 

4 

1.5 

1.92 

2.98 

4.61 

4.03 

6.53 

3.17 

l  i  able  OmUuhics) 


EXPLOSION  CONSTANT  i«c 
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Table  2  (Continued) 

Time  Measurements  from  Oscillographic  Records 


Shot 

No. 

Voltage 

(kv) 

R 

(in.) 

» , 

(/wee) 

0» 

(Msec) 

#3 

(Msec) 

04 

(/wee) 

0» 

(Msec) 

Tutor 

(Msec) 

47 

12 

12.0 

2.08 

— 

— 

_ 

— — 

... 

41 

10.5 

2.17 

3.31 

4.91 

4.25 

6.80 

3.30 

35 

9.0 

1.99 

3.50 

5.10 

444 

6.90 

3.21 

23 

6.0 

1.99 

3.21 

4.82 

4.25 

6.71 

3.02 

17 

4.5 

2.08 

3.02 

4.73 

3.97 

6.52 

3.30 

11 

3.0 

— 

— 

— 

— 

— 

— 

5 

1.5 

1.99 

3.02 

4.44 

3.97 

6.43 

3.11 

48 

10 

12.0 

1.98 

3.36 

5.00 

4.22 

6.72 

3.16 

42 

10.5 

1.98 

3.24 

4.83 

4.22 

6.72 

3.16 

36 

9.0 

2.07 

3.36 

5.00 

4.31 

6.81 

3.25 

30 

7.5 

1.99 

3.53 

5.00 

4.40 

6.64 

2.98 

24 

6.0 

1.90 

3.16 

4.66 

4.14 

6.52 

3.07 

18 

4.5 

1.98 

3.10 

4.74 

4.14 

6.50 

3.16 

12 

3.0 

2.16 

2.93 

4.48 

3.88 

6.26 

3.16 

£  (On)i 

_ 

j 

0i  =  2.07  /xsec  ft,  =  4. 15  /tsec 

0»  =  3.  12  #tsec  05  =  6.92  nsec 

$3  —  5  08  f«ee  =  3.42  /wee 


60 
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from  0.5  lb  to  55  lb  of  TNT.  Osborne  and  Taylor  (27),  using  a  millipound  of  TNT,  were  not 
able  to  report  this  phenomena,  although  they  did  show  a  departure  for  the  acoustic  law.  How- 
ever,  increasing  t,  (or  profile  broadening)  for  shock  waves  is  predicted  by  the  theory  of  Kirk¬ 
wood  and  Bethe  (12). 

A  graph  of  tmait  Figs.  12a  and  12b,  indicates  some  evidence  of  p.ofile  broadening  with 
distances.  Figure  12a  presents  tmo*  versus  distance  for  20,  18,  16,  and  14  kv,  while  Fig.  12b 
presents  data  for  12  and  10  kv.  The  strongest  evidence  for  increasing  a  time  constant  with 
propagation  is  given  in  the  20-kv  data;  the  slope  of  the  curve  draw  a  through  the  points  for 
decreasing  initial  electrical  energy  is  negative.  For  10-kv  data  it  is  very  difficult  to  establish 
that  the  slope  of  the  curve  is  positive.  Other  time  measurements,  except  05,  did  not  show  this 
be'  ,ior. 

Comparison  of  Reduced  Time  Constants  for  EWP  and  CUE  Data  —  The  time  parameter  for 
large-amplitude  CUE  events  scales,  i.e.,  forms  a  similarity  curve.  Figure  13  presents  for  CUE 
events  the  theoretical  variation  of  the  reduced  time  constant  with  the  scaled  reciprocal  distance 
in  curve  1 ,  which  is  given  by  the  expression 


jpi/a 


=  9.52 


»/* 


+  0.873 


(25) 


where  t  is  the  time  in  microseconds,  W  is  the  charge  weight  in  pounds,  and  R  is  the  gage  distance 
in  feet.  Curve  2  of  Fig.  13  is  based  on  the  experimental  work  of  Arons  (34)  and  is  given  by  the 
empirical  formula 


e 

r«/» 


=  58 


/r^V®“ 

\  R  ) 


(26) 


The  time  measurement  for  CUE  data  is  0,  in  microseconds,  and  is  the  time  required  for  the 
pressure  to  drop  from  its  maximum  value  to  (1  —  e  ’)  of  its  maximum  value  (about  67  per¬ 
cent).  The  time  measurement  t?  for  EWP  data  is  the  explosion  constant  (rr  =  2fmoj-~T,,)  ot  the 
pressute  wave  in  units  of  microseconds  time  located  at  the  maximum  of  the  curve. 


IS  -  l  .impjuvi.ti  <>(  irUiKrU  nine  cumsjnl 
v  uUtl  irtipitn.il  dnt-iiu  r  lot  (  I  I.  4!kI  IWP  il-itj 
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Fig.  14  —  Plot  of  rise  time  versus  scaled  reciprocal  distance 
corresponding  to  EWP  data  of  Fig.  13 


EWP  data  are  presented  in  Fig.  13  in  the  font,  of  points.  The  voltage  (or  energy)  associated 
with  each  point  is  identified  by  the  legend  of  this  figure.  The  EWP  data  are  obtained  from  tma x 
measurements.  Aside  from  tmax  and  0S,  the  other  0  measurements  of  Fig.  1 1  do  not  indicate  a 
negative  slope,  nor  do  they  scale.  The  EWP  data  of  Fig.  13  show  qualitative  numerical  agreement 
to  CUE  theory  and  CUE  data.  The  reduced  lime  constant  for  EWP  does  scale,  however.  An 
empirical  expression  relating  the  reduced  time  constant  to  the  scaled  reciprocal  distance  is 
obtained  from  Fig.  13. 


EWP  data  do  show  a  gentle  negative  slope;  this  suggests  profile  broadening  (increasing  tmax 
and  05  with  distance  propagation).  Upon  this  delicate  evidence,  the  nonlinear  propagation  of 
a  smali-amplilude  shock  wave  can  be  inferred,  since  profile  broadening  and  nonlinear  propa¬ 
gation  are  inseparable  constituents  of  shock-wave  propagation  (<.?.,  one  cannot  occur  without 
the  other). 

As  a  matter  of  interest,  these  results  are  contrasted  to  the  work  of  Osborne  and  Taylor, 
who  observed  a  departure  from  the  acoustic  law  but  not  profile  spreading.  In  this  study  no 
departure  from  the  acoustic  law  can  be  claimed,  but  some  evidence  of  profile  spreading  is 
observed  for  some  time  measurements. 

Figure  14  also  presents  EW'P  data  points  of  rise  time  associated  with  each  data  point  above. 


Oscillographic  Data 

Figure  15  represents  101.5  (20  kv)  joule  explosions  whose  oscillographic  records  were  taken 
at  gage  distanc  es  ranging  from  1.5  to  12.0  in.,  in  1 .5-in.  intervals.  This  set  of  records  is  represen¬ 
tative  of  the  entire  43  shots.  Close  insjrection  allows  classification  of  these  oscillograms  according 
to  their  waveform  as  type  A,  type  B,  and  type  C.  as  shown  in  Fig.  10. 

Type  .4  Oscillogram—  Tyjx*  A  oscillograms  are,  for  all  energies  (t.e.,  20,  18,  Hi,  14.  12,  and 
10  kv),  associated  with  a  gage  distance  of  1.5  in.  anti  sometimes  3.0  in.  for  higher  energy  explo¬ 
sions  and  are  characterized  bv  two  pressure  peaks.  The  first  pressure  peak  is  associated  with  the 
initial  jrortion  of  the  wire-explosion  process.  The  second  pressure  |>eak  is  approximately  one- 
quarter  the  amplitude  of  the  first,  occurs  approximately  1 7  p sec  after  the  arrival  of  the  primary 
shock  wave,  and  has  a  positive  value  of  amplitude.  The  agency  causing  this  secondary  wave  is 
not  unitpiely  determined,  but  several  jrossihle  causes  can  lx' t  tiled  out.  First,  the  secondary  wave 
cannot  be  reflected  from  the  watei  surface  alxive,  because  depths  of  1.0  to  3.0  in  lequite  a 
time  interval  of  31  to  102  /cset  ,  whic  h  is  at  variance  with  observations,  and  Ix-cause  reflections 
wo  ild  exhibit  a  negative  pressure,  which  is  not  observed  in  t \ jx-  A  oscillographs.  Secondly, 
the  possibility  of  a  bubble  oscillation  is  ruled  out  due  to  the  inertia  of  the  w.net  (production 
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Fig.  15  —  Oscillographic  records  for  1015  joule  F.WP  shots; 
time  5  p_set  per  division 


Fig.  lfi  -  (lassifu  alinn  of  pressure- 
tinu-  oscillograms 


(>l  a  secondary  pressure  wave  requires  time  of  the  order  of  milliseconds).  This  is  furthtr  ruled 
out  bv  the  experimental  work  of  Gilstein  (28),  for  example,  using  oil  instead  of  water,  which 
showed  that  bubble  formation  alone  required  several  milliseconds.  The  most  likely  contributor 
to  the  secondary  pressure  pulse  seems  to  be  a  mechanism  associated  with  the  exploding  wire 
itself;  this  would  be  the  phenomenon  of  restrike  or  reignition.  When  current  dwell  occurs 
during  the  wiie  explosion,  a  subsequent  reignition  and  expansion  of  plasma  (a  second  explosion 
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of  the  wire)  me  place.  Further  observation  of  the  secondary  shock  wave  at  increasing  values  of 
R,  due  to  either  explanation,  is  obscured  because,  while  their  amplitudes  are  different,  they 
propagate  at  approximately  the  same  speed.  To  establish  firmly  what  agency  causes  the  sec¬ 
ondary  shock  wave  is  not  possible  with  the  present  available  information;  electrical  measure¬ 
ments  of  the  explosion  process  are  necessary  for  this.  This  most  interesting  anomaly  appears, 
however,  to  be  atypical  of  CUE  events  and,  apparently,  a  property  of  EWP  events  under  the 
conditions  of  this  experiment. 


Fig.  17  —  Diagram  of  the  experiment 
showing  the  tank  dimensions 


Type  B  Oscillogram — Type  B  oscillograms  are  characterized  for  all  energies  and  for  distances 
generally  of  3.0  to  7.5  in.  b\  a  primary  shock  wave,  which  falls  f  rom  some  peak  pressure  value 
to  a  random  varying  pressure  whose  amplitude  ranges  from  0.25  to  0.15  pm,n- 

Type  C  Oscillograms -Type  C  oscillograms  are  characterized  for  all  initial  energies  and 
for  distances  ranging  from  7.5  to  12.0  in.  by  a  primary  shock  wave  which  decreases  from  some 
peak  value,  as  in  type  B  oscillograms,  but  which  differ  from  type  B  by  the  inclusion  of  a  very- 
well-defined,  steep-slope  negative  pressure  trace.  The  onset  of  mis  negative  pressure  varies 
from  10  to  23  /csec  after  the  arrival  of  the  primary  shock  wave.  The  cause  of  the  negative 
pressure  gradient  appears  to  be  the  reflection  of  the  primary  shock  wave  f'vom  the  surface  of 
the  water.  The  situation  for  EWP  in  this  study  is  described  in  Fig.  17,  where  the  wire  is  exploded 
at  point  A  and  a  shock  wave  propagates  outward  toward  point  B  and  toward  the  surface  of  the 
water  at  point  C.  At  point  C  the  wave  is  reflected  back  toward  point  B.  The  time  interval  between 
the  arrival  of  the  primary  shock  wave  at  B  and  the  onset  of  the  negative  pressure  gradient 
depends  upon  the  depth  at  which  the  explosion  takes  place.  If  Cn  is  1500  m/sec  and  d  is  2  to  3  in., 
the  time  interval  At  is  10  to  23  p sec,  respectively.  During  the  course  of  the  experiment,  because 
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no  special  effort  was  made  to  keep  the  parameter  d  constant,  type  C  oscillograms  are  observed 
for  7.5  «  R  s;  12.0  in.  for  time  intervals  of  approximately  16  ftscc.  The  tank  dimension  and 
depth  of  water  below  the  expolsion  preclude  any  other  source  for  the  observed  negative  pres¬ 
sure  gradients.  The  negative  pressure  appears  at  the  right  time  interval,  given  the  depth,  gage 
distance,  and  propagation  speed  of  the  explosion  wave.  These  parameters  are  related  by  the 
simplified  expression: 

where  At  is  the  time  interval  following  the  primary  shock  wave  of  the  reflected  shock  wave  at 
the  pressure  gage  in  microseconds,  D  is  the  gage  distance,  d  is  the  depth,  and  C0  is  taken  as  the 
speed  of  sound  for  purposes  of  approximation. 


CONCLUSIONS 

Based  upon  the  measurements  and  analysis,  several  conclusions  may  be  drawn  concerning 
the  shock-wave  parameters  pressure  and  time  and  their  comparison  to  chemical  explosives. 

1.  Peak-pressure  measurements  form  a  similarity  curve  only  if  corrections  for  switch  losses 
and  vaporization  requirements  are  taken  into  account. 

2.  As  a  result  of  the  first  conclusion,  the  principle  of  similarity  is  extended  in  terms  of  charge 
weight  six  decades  from  pounds  to  micropounds  for  the  case  of  peak  pressure,  or  similarity 
now  spans  15  decades  (i.e.,  10  9  to  10~6  lb).  A  plot  of  peak  pressure  versus  the  tube  root  of  the 
equivalent  weight  divided  by  the  gage  distance  is  w  ithin  the  same  range  of  values  given  by  Arons 
for  TNT  pressure  data. 

3.  As  a  result  of  the  first  and  second  conclusion,  the  formula  relating  the  total  initial  stored 
electrical  energy  to  its  equivalent  weight  in  TNT  is  given  by: 


W  (micropounds)  =  525 


.  /0.325  CV\?  -  Ev 


Ho 


)■ 


Here  C  is  the  capacitance  in  farads,  V0  is  potential  in  volts,  Er  the  vaporization  energy  in  joules, 
and  Hu  is  the  heat  of  detonation  of  I’NT  in  calories  per  gram.  This  formula  accounts  both  for 
the  switch  loss  of  35  percent  of  the  initial  stored  energy  and  for  the  vaporization  energy  (or 
joule  heating). 

4.  T  he  asymptotic  formula  for  peak  pressure  for  F.WP  data  is: 


Pm 


26,800 


5.  Comparison  of  the  shock-wave  energy  at  a  gage  distance  of  1.5  in.  to  the  total  initial 
electrical  energy  indicates  that,  on  the  average,  58  percent  of  the  total  initial  electrical  energy 
is  converted  to  shock-wave  energy.  This  figure  is  supported  by  independent  investigators. 

6.  Because  the  usual  techniques  for  determining  the  explosion  time  constant  in  chemical 
experiments  cannot  l>e  utilized  in  this  study,  an  approximation  based  upon  certain  assumptions 
is  made  such  that: 


TV  —  2t  iiuij-  Tii 

where  t,.  is  the  explosion  lime  constant  for  the  EWP  event,  rmai  is  the  total  rise  time  measured 
from  zero  pressure  to  peak  pressure,  and  t„  is  a  characteristic  time  constant  (the  sum  oi  the 
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amplifier  rise  time  and  the  transit  time  of  the  explosion  wave  across  the  crvstal).  By  calculation, 
Ta  is  3.41  fisec.  The  average  of  tma *  for  43  shots  is  3.42  /usee.  Therefore,  the  calculations  in¬ 
volving  Tr  and  ra  are  given  some  reason  for  confidence. 

7.  The  principle  of  similarity  (and  the  Kirkwood-Bethe  theory  consequently)  predicts 
increasing  time  duration  with  propagation.  The  time  measurements  of  the  oscillographic 
records  (8\,  0*,  03,  04,  03,  and  (max)  show  evidence  of  profile  broadening  for  0S  and  tmaxOnly. 
The  other  time  measurements  should,  but  do  not,  show  the  same  trend.  In  the  case  of  8 5  and 
tmax,  the  increase  ir,  duration  is  gentle  but  sufficiently  evident  to  be  recognized. 

8.  The  time  measurement  (tmax)  is  altered  to  give  the  reduced  time  constant  for  EWP  data 
and  is  plotted  against  its  corresponding  value  of  W'^IR,  the  scaled  reciprocal  distance.  This 
comparison  shows  that  the  reduced  time  constant  for  EWP  is  of  the  same  order  of  magnitude 
as  CUE  data,  but  that  EWP  values  are  larger.  This  comparison  shows  further  that  the  reduced 
time  constant  does  form  a  similarity  curve  as  does  CUE  data.  Finally,  because  of  the  negative 
slope,  this  comparison  shows  evidence  of  profile  broadening;  however  slight,  it  is  nonetheless 
distinguishable.  The  reduced  time  constant  is  given  by  the  empirical  formula: 

(■ tM 1/3 )  =  70  (W'l'lR)-0*1 

9.  The  data  taken  at  R  =  1.5  in.  and  some  at  R  =  3.0  in.  show  a  secondary  peak  pressure 
having  a  positive  amplitude  of  0.30  pm  of  the  first  peak-pressure  amplitude.  A  speculative  guess 
as  to  the  origin  leads  the  author  to  believe  that  currert  dwell  and  subsequent  reignition  (a  second, 
but  lesser  explosion)  take  place.  Correlated  electrical  and  optical  measurements  (not  made  in 
this  study)  would  resolve  this  question. 

10.  Much  of  the  data  taken  at  R  3*  7.5  in.  show  a  large  negative  pressure  (a  reflected  shock 
wave)  incident  roughly  17  /csec  after  the  peak  pressure  of  the  explosion  wave.  Since  the  tank 
dimensions  preclude  reflections  from  the  sides  and  bottom  of  the  tank,  this  reflection  must  be 
from  the  air-water  interface.  This  conclusion  is  supported  by  calculation. 
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